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1. Introduction 

1.1. New laws from old ones 

If general relativity is "probably the most beautiful of all existing physical theories" 
|14| . it is certainly thanks of its geometric character, which reduces the dynamics 
of test bodies in a gravitational field to pure kinematics^ This feature makes it an 
unprecedented heuristic machine: to uncover the effect of gravity on a given physical 
phenomenon, consider the old law which describes it in the absence gravity, phrase it 
in kinematical terms (times lapses and distance intervals), replace "time" by "proper 
time" and "distance" by "proper distance" , and read off the new law in the presence 
of gravity. For instance, the Fermat principle states that, in the absence of gravity, 
light rays follow the paths which extremize time. Then general relativity immediately 
generates a new law from Fermat's principle: in the presence of gravity, light rays 
follow the paths which extremize proper time. The bending of light in the presence of 
spacetime curvature follows immediately from this new law. 

It is interesting to note that the generative character of general relativity is 
actually more general than general relativity itself. Indeed, it relies neither on local 
Lorentz invariance, nor on the absence of a preferred foliation of spacetime, nor on 
the Einstein equation, nor even on the relationship between stress-energy density and 
spacetime curvature. This is particularly clear in the example mentioned above: the 



I "Kinematics" has several, inconsistent, meanings in the physics literature. Here, we use this term 
to intend the description of a phenomenon purely in terms of space and time. 
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spacetime curvature responsible for the bending of light could come from a refractive 
index gradient as well as from the vicinity of a massive star. In other words, whether 
spacetime curvature is fundamental or effective is irrelevant to the generation of new 
law from old ones: all that matters is that a curved spacetime, seen through geodesic 
coordinates, appears flat. 

This circumstance is all the more important that a growing number of condensed- 
matter systems are now understood to behave as effective, or analogue, spacetimes. 
Besides Gordon's refractive optical medium, one can mention Unruh's dumb hole (a 
supersonic fluid flow), but also corrugated graphene sheets, Bose-Einstein condensates, 
slow light systems, superfluids, metamaterials, etc. (See [1] for an updated review 
of analogue gravity.) The fruitfulness of this connection between condensed- matter 
physics and general relativity goes both ways: gravitational analogues provide valuable 
model systems to emulate otherwise out-of- reach relativistic phenomena |25) : vice 
versa, the geometric setup of general relativity sheds a new light on venerable fields 
such as optics |16 ] , hydrodynamics [121 - or diffusion phenomena, as this paper intends 
to demonstrate. 

1.2. The Tolman-Ehrenfest law of thermal equilibrium 

That the "new law from old ones" principle does not restrict to the realm of mechanics, 
or electrodynamics, is demonstrated by the early history of general relativity. Indeed, 
one the first problems which Einstein analyzed in terms of gravitational redshift 
was one of thermodynamics: the problem of finding the equilibrium temperature 
distribution T* in a static gravitational field. As early as 1912, that is three years 
before the completion of general relativity, he speculated that, because in a curved 
spacetime proper time does not run at the same rate in different places, T* should not 
be homogeneous [9^. This remarkable intuition was put on firm ground by Tolman 
and Ehrenfest, who showed that 

xT* = const., (1) 

where x = ("C^Ca)^^^ is the redshift factor and a timelike Killing vector f2T. 

From our perspective, the Tolman-Ehrenfest relation ([T]) is a prototype of these 
kinematical laws which follow directly from the geometrical setup of general relativity. 
In their original paper, however, Tolman and Ehrenfest gave it a complicated 
dynamical proof, which relied both on the Einstein equation and the equation of state 
of thermal radiation. Several authors later pointed out this anomaly, and proposed 
more minimalist derivations [6l[23]. One, due to Rovelli and the author, goes like this 

m- 

• In the non-relativistic canonical ensemble, the equilibrium temperature T* can 
be computed as the rate of the modular flow generated by a thermal state p (the 
Hamiltonian flow of Inp) with respect to time. 

• Hence in a stationary spacetime, T* can be computed locally as the rate of the 
modular flow generated by a thermal state p (the Hamiltonian flow of In p) with 
respect to proper time. 

• By stationarity, the modular flow generated by p is proportional to the Killing 
flow, hence T*{a) cx dt{cF) / ds{a) , where t{a) is the Killing parameter at a spatial 
point cr and s{a) the local proper time. 

• By definition, this ratio is the inverse of the redshift factor x(c)i hence T*(cr) oc 
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In other words, the equihbrium teniperature distribution T* is such that the 
combination xT* satisfies the non-relativistic equihbrium criterion. Can this reasoning 
be extended to non-equihbrium processes such as heat conduction, or more general 
diffusion processes? How far does the "new laws from old ones" principle lead us 
when we leave the equilibrium regime? 

1.3. The stochastic route to diffusion 

As already emphasized, the key step to address this question is to frame the diffusion 
problem in kinematical terms. Unsurprisingly, the most direct route to that effect is 
the one laid down by Einstein himself, in his 1905 work on Brownian motion 8J: write 
the master equation for a stochastic process, and derive the corresponding Fokker- 
Planck equation in the diffusive limit. 

In this stochastic approach, the basic concept is that of transition rates: the 
probability that a random walker will jump from one position to another per unit time. 
These are kinematical in nature (they are expressed in terms of distance and time), 
and are therefore easily amenable to the "new laws from old ones" principle. Precisely 
by this token, we will derive in this article the curved-spacetime generalization of 
the Fokker-Planck equation, applicable to any kind of diffusive transport (atomic and 
molecular diffusion, photon diffusion, thermal and electronic conduction, etc.), in any 
kind of fundamental or analogue spacetime. In the case of pure Brownian motion and 
static spacetimes, it reads 

eVa(xp) = KA(xp), (2) 

where p is the probability density of Brownian motion, = /x the hydrostatic 
4-velocity, A the spatial Laplace-Beltrami operator and n the diffusivity. Note that 
equation ^ is nothing but the standard diffusion equation, with p replaced by xPi as 
in the Tolman-Ehrenfest relation - plain and simple. 

Although (to our knowledge) it was never written in this form, equation ^ 
is actually well-known in relativistic hydrodynamics: it is the phenomenological heat 
equation of Eckart [7] and Landau and Lifschitz [15] . Here, it is derived from stochastic 
mechanics, rather than postulated to satisfy the second law of thermodynamics - just 
like Einstein derived the diffusion equation postulated by Fourier and Fick. 

H. On causality 

The notion of "relativistic Brownian motion" has been discussed by many authors in 
the past decades [5]. What is usually meant by this expression is a formulation of 
Brownian motion that is consistent with special-relativistic causality. Our goal here 
is different: we wish to understand the effect of a non-trivial spacetime geometry on a 
diffusion process. 

As is already apparent from our equation ([2]) , which is parabolic and hence permits 
superluminal propagation, we do not attempt to include causality in our framework. 
Instead, our setup is the following. We consider a fluid flow in spacetime - the bath 
within which the stochastic process takes place - and fix the associated orthogonal 
foliation of spacetime. In this foliation, each hypersurface is the "instantaneous space" 
relative to fiducial observers comoving with the fiow; at each instant i, a random 
walker dragged by the fiow can jump from one point to another point of the same 
hypersurface. 
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Thus, the sole difference between our approach and the non-relativistic theory is 
the inclusion of a non-trivial spacetime geometry, viz. the presence of a gravitational 
field. That is to say that we consider the c — )■ oo limit of a fully relativistic theory 
of Brownian motion, in the same sense in which Newton-Cartan gravity |101 [5] is the 
c ^ 00 limit of general relativity. 

Of course, this acausal approach is at variance with the textbook notion that 
special relativity precedes general relativity. Yet, it is by no means unreasonable. It 
simply expresses a physical approximation, namely that the relaxation of the bath 
occurs on time scales much shorter than that of the diffusion process itself. This 
regime is well-known in the context of kinetic theory, where it is referred to as the 
"hydrodynamic limit" , and is also the regime in which Eckart's dissipative relativistic 
hydrodynamics appliesl§| We refer the reader to [3] for an interesting discussion on 
why an acausal equation like ^ is not inconsistent with microscopic causality. 

1.5. Results 

Besides deriving the master and Fokker-Planck equations for stochastic processes in 
curved spacetimes, our results in this paper are the following. 

• We provide a microscopic justification to Eckart's heat-flux ansatz, and extend it 
to more general diffusion processes. 

• We generalize the Tolman-Ehrenfest relation to non-equilibrium stationary states, 
with arbitrary boundary conditions. 

• We compute the gravitational corrections to the mean squared displacement of 
Brownian motion in static isotropic spacetimes. 

The last item is particularly interesting. In a curved spacetime, the usual scaling law 
(x^)t oc t holds only in the t asymptotic limit. At later times, spacetime curvature 
corrections show up and modify the growth rate of {x'^)t- This suggests that diffusive 
transport in gravitational analogues could perhaps be tailored, by tuning the effective 
metric coefficients [221 . 




1.6. Plan of the paper 

The paper is organized as follows. Sec. [2] consists of preliminaries on the 
D+1 formalism for relativistic hydrodynamics and on the non-relativistic theory 
of stochastic processes. Our theory of stochastic processes in curved spacetimes is 
developed in sec. [31 and the limit case of Brownian motion is studied in sec. [3| In sec. 
[SJ we obtain a small-time asymptotic expansion for the mean squared displacement of 
Brownian motion in static isotropic spacetimes. Our conclusion follows in sec. [6| 

2. Preliminaries 

Throughout this paper, we consider a. {D + 1) dimensional spacetime with signature 
(— -I- -I----). (We keep D unspecified to include lower dimensional analogue 

§ Eckart's theory is often considered "unacceptable" because of its acausal character, and its alleged 
instability [11] . As far as the author can see, this judgement is completely misled; Eckart's heat 
equation is a perfectly well-behaved parabolic PDE, whose status with respect to a fully relativistic 
dissipative hydrodynamics is the same as that of Newton-Cartan gravity with respect to general 
relativity: an excellent approximation in most physical situations. See |13) for a mathematical 
argument to this effect. 
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spacetimes in the discussion.) We denote V the spacetime Levi-Civita connection, 
and a, 6, c..., i, j, ... are abstract indices. 

The standard references for general relativity and the D + 1 formalism are [T71I28) : 
stochastic processes and Fokker-Planck equations are exposed in [371 HO] ■ 



2.1. The D+1 formalism 

Consider a relativistic fluid with velocity Assume its flow is irrotational, viz. 

U[aVbMc] = 0. (3) 
Then, according to the Frobenius theorem, there is a foliation of spacetime by 
hypersurfaces Ef orthogonal to Furthermore, the slices St are the level sets of 
a time functions i : ill — > M such that 

Ua = -NVat (4) 

for some non-negative function N . The function TV is called the lapse function^ and 
the slices Et have the interpretation of "instantaneous space" relative to observers 
comoving with the fluid. In the following, we will denote a a flow line of u° (a "spatial 
point"), and at its intersection with Ei. 

The intrinsic geometry of the spatial hypersurfaces Et is coded by the induced 
metric 

hab = gab + UaUb, (5) 

and its associated covariant derivativ^ Da and Laplace-Beltrami operator A, while 
their embedding in spacetime is measured by the (symmetric) extrinsic curvature 
tensor 

Kab = Valifc. (6) 

The trace 9 — — Vqu" of the extrinsic curvature tensor is called the expansion 
scalar. It measures the fractional rate of change of an infinitesimal volume 5V about 
a spatial point along the flow, viz. 

The factor above converts the proper time along the flow to the global time 
coordinate t. 

A situation of particular interest is the hydrostatic equilibrium: the vector 
ya^ „ —^a^jy jg ^jj^gjj Killing, 1.6. generates timelike isometrics. In this context, 
the lapse function N is usually denoted x, and called the redshift factor. It satisfles 
■""VaX = 0, and gives the acceleration = u'^VcU^ of the flow by 

a'^V^lnx- (8) 

Moreover, the time-time component of the Ricci tensor E = RabU°'u^ (sometimes 
called the Raychaudhuri scalar), is given in this case by 

E = Dba^ + aba^. (9) 

In general relativity, this scalar is tightly related to the local mass-energy density, by 
virtue of the Einstein equation. We will see that E plays an interesting ro in diffusion 
phenomena. 

II The covariant derivative Da associated to h^b acts on a tensor field T"^ t"...;, according to 
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2.2. Markov processes 

Let E be a Riemannian manifold with metric hij and covariant derivative Di, 

representing a curved space, and denote (Tf G E the instantaneous position of a random 
walker at time t. In the Markovian setup, we assume that at completely determines 
its later positions at' {f > t), according to transition rates r{a — > a'). By definition, 
these arc such that the elementary probability for the walker to jump from a volume 
dV{a) about a G E to a volume dV{a') about a' G E in time dt is given by 

T{a a')dV{a)dV{a')dt. (10) 

As a rule, the transition rates are implicit functions of the metric hij. 

Let Pt{a) denote the probability density that the walker is in neighborhood of a 
at time t, i.e. at = a, and 

jtia^a')=pt{a)r{a^a'). (11) 

the corresponding probability fluxes. Balancing the incoming and outgoing fluxes at 
a, we can immediately write the evolution equation for pt as 

dtPticr) = dV{a') (it(a' ^ a) - jt{a ^ a')) , (12) 

i.e. 

dtPticr) = dV{a') {pt{a')T{a' ^ a) - Pt{a)T{a ^ a')) , (13) 

where dV{at) is the Riemannian volume element on E^. This integro-differential 
equation is known master equation, and the operator A4 such that dtPt = M.pt as 
the master operator. 

In this stochastic framework, the notion of equilibrium state has a clear-cut 
definition: a steady-state solution p* is an equilibrium distribution if the corresponding 
probability fluxes cancel pairwise, i.e. 

p*ia)T{a ^ a') =p*ia')r{a' ^ a) (14) 

This condition is known as the detailed balance condition. 

Under certain regularity conditions for the rates F, one can show that the paths 
{at) are discontinuous: for this reason one often speaks of jump processes in this case. 
The situation changes in the limit where the jumps become inflnitely frequent and 
short-ranged (with respect to some relevant coarse-graining scale). Then F becomes 
distributional, and the master operator M reduces to its second-order truncation L in 
a moment expansion, reading 

Cpt = -D,{w\pt) + ]^D,D,{w^pt). (15) 

Here w\ is a vector field on E, the drift vector, and W2 a symmetric and positive- 
definite rank-2 tensor field, the diffusion tensor. Note that the transition rates F are 
related to C according to 

T{a' ^a)=C5{a',a), (16) 

where 5 is the Dirac distribution on E and C acts on the a' variable. Stochastic 
processes described by such Fokker-Planck equations are called diffusion processes. 

The simplest example of such a diffusion process is Brownian motion, for which 
(by definition) = and W2 = 2Kh^^ for some positive constant k. The 
corresponding Fokker-Planck equation dtPt = iCpt is the well-known diffusion equation 

dtp = KApt. (17) 
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3. Master and Fokker-Planck equations in curved spacetimes 

In this section we describe the curved spacetime generahzation of the master and 
Fokker-Planck equations for Markov processes. 



3.1. Markovian setup 

Consider a Markov process defined by stationary transition rates T{a' — >■ a), depending 
parametrically on a Riemannian metric hab- In the case of Brownian motion, for 
instance, r{a' a) = kAS{(t,(t'), with A the Laplacc-Bcltrami operator associated 

to hab. 

Following the general "new law from old ones" ansatz, we take this process as 
defining the instantaneous dynamics of a random walker in spacetime, in proper time. 
In other words, given an irrotational flow u", we consider the associated orthogonal 
foliation (St), evaluate F on the induced metric /ia60 and assume that the probability 
that a random walker carried by the flow u° will jump from the position at to the 
position cTj in proper time ds{at) is given that 

r{at a't)dViat)dVia't)dsiat), (18) 

where ds{at) is the proper time along a. 



3.2. Master equation 

Now, to write the corresponding probability equation, which is necessarily global, we 
must convert the proper time s((Tt) in (|18p into the time coordinate t. This is achieved 
thanks to the lapse function N , as 

ds{(Tt) = N{at)dt. (19) 

Hence, we can rewrite (|18l) as 

T{at a[)dViat)dVia[)Niat)dt. (20) 

Denoting p{at) the probability density of the stochastic process, the probability flux 
is therefore 

jiat ^ a[) = N{at)p{<Jt)T{at ^ a'^). (21) 

This expression is physically intuitive: where proper time runs faster (high N), the 
walker jumps more frequently (high j). 

From this simple argument, we get that, if A4 is the master operator associated 
to the rates F, the right-hand side of the curved-spacetime master equation should be 
M{Np), i.e. 

£ dV{a',) (^N{<j[)pt{<j[)r{a', ^ at) - N{at)pt{a)r{at ^ a[)) . (22) 

A moment of reflection shows that the left-hand side of the master equation should 
also be modified in a curved spacetime. Indeed, recall that in a curved spacetime, the 
time- variation of an integrated density does not coincide with the integral of the time- 
derivative of the density: it Vj is a region in Ef, then 

^ / dV{at)pt{<Jt) + [ dV{at)dtPt{<7t). (23) 



^ If spacetime is not static, this makes the transition rates implicit functions of time. 
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This is due to the fact that the vohime element dV{(Jt) itself depends on time. The 
correct formula follows from the relationship ([7]) defining the expansion scalar, and 
reads 

j^dV{at)pt{at)^ j^dV{at)(dtPt{at) + NBpt). (24) 

Shrinking the volume Vt down to zero, we thus find that the left-hand side of the 
master equation should be dtP + NOp instead of dtPt ■ 

Combining both insights, we find that the master equation in a curved spacetime 
with lapse N and expansion 6 is 

dtp + NOp^ M{Np). (25) 

It is easy to check that this equation conserves the total probability Jj, dV{(Jt)pt{(Jt), 
as it should. 



3.3. Detailed balance condition 

Note that, in the case of static spacetimes {0 = and TV = x is the redshift factor), we 
can read off from (I^TI) the generalized detailed balance condition: for an equilibrium 
distribution p* , the probability fiuxes cancel pairwise if 

r(a' ^ a)x{a')p*ia') = r{a -> a')x{a)p* (a). (26) 

Hence, the product XP* must satisfy the usual detailed balance condition defined by 
the rates r(CT — ?► a'), instead of p* itself, as in the non-relativistic case. This is the 
stochastic counterpart of the Tolman-Ehrenfest relation ([1]), where xT* itself satisfies 
the usual homogeneity condition instead of T* . 



3.4. Diffusive limit 

Assume from now on that the stochastic process is of diffusive type (or can be 
approximated by on43) and denote C the Fokker-Planck operator defined by the rates 
r, as in ([15]). Then from ([25]) it follows immediately that the Fokker-Planck equation 
reads 

dtP + Nep = L{Np), (27) 

i.e. 

dtp + NOp = -Da{wlNp) + ]^DaDk{wfNp) (28) 

where and W2'' are the drift vector and diffusion tensor associated to the rates F, 
as in sec. 12.21 This is the curved- spacetime Fokker-Planck equation. 

Note that (|28p can be given a more hydrodynamical flavor, by replacing the 
unphysical derivative dt by the convective derivative w^Va, which evolves the 
probability distribution in proper time rather than in coordinate time; it then becomes 

u^v^P + ep = + iDM<m^ (29) 

al^-r p N 2 N ^ ' 

This equation is the main result of this paper. 



+ We recommend van Kampen's note [26l for a discussion of the applicability of this approximation. 
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In this section we focus on the properties of Brownian motion in curved spacetimes. 

The general-relativistic diffusion equation 

We saw in sec. [Hthat Brownian motion is characterized among diffusion processes by 
the vanishing of the drift vector, = 0: ^-^d by = k/i"'', with k the diffusivity. 
The corresponding Fokker-Planck equation is therefore 

dtp + NOp = kA(7Vp) (30) 

or 

U-VaP + ep=K^^. (31) 

The remainder of this article is concerned with the properties of this curved- spacetime 
diffusion equation. 

4-2. Comments in the hydrostatic case 

Consider the hydrostatic case, where ([5T|) reduces to 

u^VaP = K^^. (32) 
X 

Several comments can be made about this equation. First, since u^VaX = Oi this 
equation indeed coincides with ([2]), as announced in the introduction. Second, using 
the relation Of, = Df, log x between the acceleration of the congruence ab and the 
spatial gradient of the redshift factor, the equation (|32|) can be reorganized as 

(^ - 2Ka')VbP = kAp + kEp, (33) 

where E is the Raychaudhuri scalar. In addition to the usual diffusion term Ap, this 
equation contains two remarkable terms, which have no analogue in the non-relativistic 
diffusion equation. 

• Drift. The term 2Ka''Vt,p is a drift term. Unlike the drift term in the classical 
Fokker-Planck equation ((T5|) . it vanishes in the limit k ^> 0, and is therefore a 
genuine effect of diffusion. 

• Source. The term nEp, where E ~ Di,a^ + Oba^, is a source term. It implies that 
the probability density appears to comoving observers as sourced by (k times) 
the Raychaudhuri scalar E^ 

Both terms, which result from the non-homogeneity of x in space, can be interpreted 
as stochastic gravitational redshift effects. 

* That is not to say that the total probabihty is not conserved; we saw that it is. 



Diffusion in Curved Spacetimes 



10 



4.3. Derivation of Eckart's constitutive relation 

Another interesting consequence of our stochastic derivation of the diffusion equation 
in curved spacetimes is the vindication of Eckart's phenomenological constitutive 
relation for the heat flux in general relativity f^: 

= ~k{D^T + Ta^). (34) 

This relation was postulated by Eckart on the basis of thermodynamical arguments, 
and can be used to write the relativistic heat equation as 

u-'VaT +{Db + ab)q^ = Q. (35) 

Consider now the diffusion equation (1321) for the probability density of Brownian 
motion in a static spacetime, and compare it with Eckart's heat equation psp : they are 
the same. In other words, we have reappraised Eckart's heat equation as a probabilistic 
equation - just like Einstein did with Fick's diffusion equation. 

^.-i- The non- equilibrium Tolman-Ehrenfest condition 

Another straightforward consequence of the lapsed heat equation is the generalization 
of the Tolman-Ehrenfest condition to non-equilibrium steady-state solutions. Indeed, 
we see from ([5(Il) - PT|l that the the steady-state solution r°° is given by ([Ij only in 
the absence of an external forcing on the boundary; in general, it satisfies instead 

A(xr°°) = 0. (36) 

Hence, steady-state solutions can be described as T°° — ^/Xi where is a harmonic 
function. (Equilibrium distributions corresponding to the case V' = const.) To our 
knowledge, this characterization of steady-state temperature distributions in static 
spacetimes was not derived before. 

5. Corrections to the mean square displacement 

In this section we compute the gravitational corrections to the mean square 
displacement of Brownian motion as a function of time. 

5.1. Assumptions 

To avoid dealing with the drift effect mentioned in sec. 14. 2[ we assume from now 
on that space is radially symmetric about o, the origin of the Brownian motion. We 
also assume that the metric is quenched, i.e. evolves at a much slower rate than 
the diffusion process itself. In this approximation, the lapse function TV and spatial 
geometry hab are essentially independent of t, and the expansion scalar 9 is negligible 
with respect to the (inverse) diffusion time, hence ([50]) reduces to 

dtKt = ^HxKt). (37) 

Hereafter, we shall denote S the time-independent spatial section, and (T, the 
pairing between a distribution T and a test function on S. We also assume (without 
loss of generality) that N{o) = 1. Finally, we disregard the possible existence of cut 
loci in S, and effectively restrict our attention to a convex normal neighborhood of o, 
where the (spatial) Riemannian distance p((t) = d{a, 0) is a smooth function of a. 
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5.2. Green function and mean square displacement 

The most significant observable of Brownian motion is tlie mean square displacement 
(MSD). It is defined as the expected value of the squared distance between the position 
of the Brownian walker at time t and its initial position: 

{p'h^iKup'). (38) 

Here Kt is the Green function (or heat kernel) of the diffusion equation ([57)) . namely 
the solution with initial condition 

\im Ktia) = d{a,o), (39) 

where 5{a, a) is the Dirac distribution on the spatial slice E with support at o. (Note 
that, with the definition (I55|) . the MSD is measured as a function of the t coordinate, 
which coincides with proper time only at the origin o: unlike the non-relativistic 
situation, there is no global physical time parameter in a curved spacetime.) 



5.3. Asymptotic expansion of the MSD 

Let us denote V the differential operator KAg(x-). Then the equation (|57)) can be 
solved formally as 

Kt{a) = e'^S{a, o) = V -:V^S{a, a). (40) 
— ' n\ 

The MSD, in turn, can be computed by evaluating this distribution the squared 
distance function p^. To this effect, note that 

{V6,p^)^{S,V^p')=V^p'{o), (41) 

where = Kx^q is the formal adjoint of V. Hence 

(p2), = (if„p2) = f;^(pt)V(o), (42) 



71! 

n=0 



I.e. 



{P% = ^-^ixArp'io). (43) 

This formula provides the asymptotic expansion of the MSD in the small time limit 
t — > 0. Up to second order in t, it gives 

ip^)t^2.Dt{l+{^-^yt + Oit^)}. (44) 

To arrive at this expression we used the geometric identities Ap'^(o) = 2D and 
A2p2(o) = _4i?(^)(o)/3ll At this order, we thus see that diffusion is enhanced by a 
convex lapse profile about o, and by negative spatial curvature. 

tl The higher order terms involve higher derivatives of the squared distance function, which can also 
be expressed in terms of local curvature invariants (4] I18j . 
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5.4- The backward equation 

Note that the expansion (|^^ can be resummed formally as 

{p% = e'^'p'io). (45) 

Thus, the MSD {p^)t can also be obtained as the solution ut to the adjoint, or backward, 
equation 

dtut = KxAut (46) 

with initial condition wo(o') = P^(o')- This differential formulation can be useful to 
obtain the MSD in concrete situations, by means of a numerical integration of (|46p. 



5.5. Two examples 

We close this section with two explicit examples where the MSD is altered by the 
spacetime geometry. The first one is the simplest general-relativistic star model, and 
the second one is inspired by condensed-matter gravitational analogues such as graded- 
index optical fibers. 

• Schwartzschild's constant density star. This is a static solution of the Einstein 
equation with uniform mass-energy density. It has two parameters R and M, the 
radius and mass of the star respectively. (See |19j for the explicit expression 
of the line element.) If o is the center of the star (r = 0), one computes 
Ax(o) = 3GM/R^ and R^{o) = 12GM/R^, and therefore 

{p')t^6nt(l + ^Kt + 0{t')). (47) 



6i?3 

Thus, a Brownian motion initialized at the center of the star spreads slightly 
faster than in flat spacetime. This result might seem paradoxical: doesn't gravity 
attract! Recall however that Brownian motion takes place within the stellar 
medium, which is not free-falling but static. The infinitely frequent collisions 
between this medium and the Brownian particle prevent the latter from falling 
down to the center of the star - on the contrary, we see here that they actually 
increase the MSD. Furthermore, this effect is small: a simple computation shows 
that the Brownian motion hits the surface of the star ((p^)^^ ~ R) long before 
the corrective term {GM/6R^)Kt becomes of order 1. 

Parabolic lapse profile. Interestingly, this speed-up effect can be emulated, and 
amplified, in a gravitational analogue with fiat spatial geometry and a parabolic 
lapse profile 

Xip) = 1 + ep^R^ (48) 

Here e = ±1 indexes the convexity/concavity of the profile. Such lapse profiles 
arise e.g. in graded-index optical fibers, or in Kerr media controlled by intense 
laser pulses. (In these optical contexts, the lapse function is nothing but the 
inverse of the refractive index.) Moreover, this case has the advantage that the 
asymptotic expansion (1421) can be resummed explicitly. Indeed, we have 

Ax = eAp^/R"^ = 2eD/R^, (49) 

hence the formula (H^ gives 



°° ^.Uf-n ,o 7-). „_i 



=1 
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i.e. 

In the convex case (e — 1), the MSD therefore grows exponentially with time, 
while in the concave case (e = —1), it sfows down and eventually reaches the 
finite limit on the "infinite redshift surface" {x = 0}. This effect becomes 
significant on the time scale R^/k. Materials where R can be tuned at will could 
therefore provide experimental benchmarks for the results discussed in this paper. 

6. Conclusion 

Let us conclude. From a theoretical standpoint, our reasoning in this paper is very 
straightforward: it simply consists in incorporating gravitational redshift and spatial 
curvature effects into the standard master equation for a Markov process - in short, 
Einstein 1905 [S] amended by Einstein 1912 [9^. 

Simple as it is, though, this approach has allowed us to derive Eckart's constitutive 
relation for heat transfer, to generalize it to non-thermal diffusion processes, and to 
compute the gravitational correction to the diffusion square-root law. In particular, 
we have obtained the general asymptotic expansion of the mean-squared displacement 
in static spacetimes, and concluded from two explicit examples that experiments are 
more likely to reveal such corrections in analogue gravitational systems. Given the 
ubiquity of diffusion phenomena in condensed-matter physics, we are hopeful that 
these results will prove useful in applications. This would confirm - if that was needed 
- that general relativity remains as fertile as ever. 
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